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Abstract 

We show that excitations of physical interest of the heavenly equa- 



tion are generated by symmetry operators which yields two reduced 
' equations with different characteristics. One equation is of the Liou- 

' ville type and gives rise to gravitational instantons, including those 

found by Eguchi-Hanson and Gibbons-Hawking. The second equa- 
tion appears for the first time in the theory of heavenly spaces and 
■ provides meron-like configurations endowed with a fractional topolog- 

. ical charge. A link is also established between the heavenly equation 

I and the socalled Schroder equation, which plays a crucial role in the 

'nJ" ' bootstrap model and in the renormalization theory. 

a^ . PACS 04.20J, 02.30 



An important statement concerning self-dual, euclidean, Einstein spaces 
Q^l with one rotational Killing vector (heavens), is that these are completely de- 

termined by a real scalar field obeying a nonlinear partial differential equation 
which formally coincides with a continuous version of the Toda lattice [||]. 
This equation can be written as 



^ , l^xx ~l~ ^yy ^(e "^zz (l) 

where the signature factor k = ±1 has been introduced for later convenience. 

Equation (P appears in a variety of physical areas, ranging from the 
theory of Hamiltonian systems to general relativity. In the last context, the 
metric tensor related to Eq. (|l|) is 



{e^iidx)' + idyf] + {dzf} (2) 



1 2 

+ — [e{—Uydx + Uxdy) + dr] 
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where e = ±1 for self-dual and anti-self-dual cases, respectively. 

Heavenly spaces have been studied by several authors in different frame- 
works 101 . The main aspects of these investigations concern i) the problem 
of the classification of the Riemann metrics admitting at least a Killing vec- 
tor field; ii) the discover of extended conformal symmetries, i.e. the large n 
limit of the Zamolodchikov Wn algebra [|l|; Hi) the search of exact solutions 
of physical significance, such as self-dual gravitational instantons [Q, whose 
role is important in quantum gravity 0]. 

In this note, we focus the attention on both the symmetry structure 
and the determination of new exact solutions of Eq.(l). In regard to the 
last point, even particular solutions are hard to find. An effective method 
to shed light on the symmetry properties of Eq.(l) and, correspondingly, 
to construct explicit configurations is based on the reduction approach 
which exploits group-theoretical techniques. Following these ideas, here we 
show that Eq.(l) allows an infinite-dimensional Lie algebra, a realization of 
which is given by generators whose coefficients are harmonic functions in the 
variables x, y. Some generators, say X„, obey the Virasoro algebra £ without 
central charge. The remaining generators, say X„, which are constructed by 
the harmonic functions conjugate to those related to the Virasoro operators, 
satisfy a set of commutation relations involving the presence of both X„ and 

Furthermore, we find two reduced equations which furnish two classes of 
solutions to Eq.(l) with quite different properties. One of them, Eq. (|l^) , pro- 
vides meron-like gravitational configurations which allow a fractional topo- 
logical charge. The second equation, ( p7b|) , yields instanton configurations 
(with integer topological numbers) containing the Eguchi-Hanson and the 
Gibbons-Hawking gravitational instantons as particular cases. Finally, an 
unexpected connection is shown between the heavenly equation and the so- 
called Schroder equation which is important in the statistical bootstrap model 
and in renormalization theory [Q. 

It turns out also that after a suitable change of basis, the elements 
Xq, Xi, Xq, Xi of £ fulfil the commutation relations of the al- 
gebra of the proper Lorentz group S0{3, 1), which is locally isomorphic to 
SL{2, C) 0. Thus, the appearance of SL{2, C) in the theory of the heavenly 
spaces is a direct consequence of the symmetry structure of Eq.(l). This 
property is also partaken by the Newman-Penrose equations, i.e. the null 
tetrad formulation of the Einstein field equations which are invariants under 
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the group SL{2,C). Therefore, the choice of SL{2,C) as the gauge group 
used in the local gauge approach to Einstein's general theory of relativity 
arises in a natural way by symmetry considerations. 

The symmetry algebra of Eq.(l) can be obtained applying the procedure 
of Ref.0. The generator of the corresponding symmetry group has been 
found in part with the help of a computer, by using the symbolic language 
REDUCE 1^. Our result is that Eq.(l) admits an infinite dimensional Lie 
group of symmetries, that is a local group G of transformations acting on the 
independent variables (x, y, z) and the dependent variable m, with the prop- 
erty that whenever z) is a solution of Eq.(l), then u' = {gou){x', y', z') 
is also a solution for any g & G. The symmetry algebra underlying G is rep- 
resented by the vector fields 

X{f) = fd, + fdy- 2/,9„, (3a) 

Y = zd, + 2du, (3b) 

T = dy, (3c) 

Z = d,, (3d) 

where / = f{x,y) is an arbitrary harmonic function, and / = f{x,y) is its 
conjugate. 

We observe that the translation operator along the axis x is contained in 
(|3a]) , in the sense that dx = -^(1)- The generators (^ obey the commutation 



relations 

[Xif),Xig)] = Xifgx - gfx + gh - fg.), (4a) 

[X{f),T] = -X{fy), [Z,Y] = Z, (4b) 

[T,Z] = [X{f),Z] = [T,Y] = [X{f),Y] = 0. (4c) 

The presence of (harmonic) arbitrary functions tells us that the symmetry 
algebra (|^) is infinite-dimensional. This algebra can be embedded into a 
loop algebra of the Virasoro type without central charge. This can be shown 
by assuming that / and g are of the form = — e"^ cos ny and gm = 
— e™^ cos my {n,m E Z). In such a way, we build up a sub-algebra of the 
symmetry algebra (^) with the basis Y,T, Z, and 

Xn = X{fn) = e"^ (— cosnydx — sinnydy + 2n cos nydu) , (5a) 

Xn = X{fn) = e"^ (— sin nydx + cos nydy + 2n sin nydu) ■ (5b) 



2 



The commutation relations 



[X„, Xra] = {n- m)Xn+m, [Xn, T] = nX„, (6a) 

[Xn, Xm] = {n- m)Xn+m, (6b) 

[Xn, Xm] = -{n - m)Xn+m, (6c) 

[Xrr, T] = -nX„, Y] = Z] = 0, (6d) 



hold, together with (^) and (^). 

The existence of the Virasoro-type algebra (^ reflects the invariance prop- 
erty of Eq.(l) under a conformal transformation. In other words, if u{x, y, z) 
is a solution of Eq.(l), the function u = u{x',y',z) given by 

u = u{x,y,z) +\n{U^^, + V^,), (7) 

is also a solution, where x = U{x',y') and y = V{x',y') are any pair of 
conjugate harmonic functions. 

Equation (|^) is a direct consequence of the invariance of Eq.(l), written 
in the complex form 

Ad^df^u = K^^e", (8) 

under the transformations 

V = f{0, V = m, u = u-\n{f'f'), 

where r] = x + iy, f] = x - iy, = ^{d^ - idy), df^ = \{d^ + idy), f = /^(C), 
/' = /^(C) and f{() is an arbitrary holomorphic function of ( = x' + iy'. 

It is noteworthy that the commutation relations (^ contain those of the 
proper Lorentz group 5*0(3, 1). This can be seen by putting 

Ji = ^(Xi-X_i), J2 = -^(Xi+X_i), 





), J2 = - 


^3 


= -iXo, 




), P2 = - 


^3 


= -iXo. 



Pi = ^(Xi-X_i), P, = -i(X, 
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Then, 



[i/j, t/j] (^ijkJki (^ijkJki 

[Ji^Pj] =^ijkPk- (9) 

Some important solutions which can be identified with meron-hke config- 
urations and gravitational instantons emerge applying the method of symme- 
try reduction to Eq.(|I]) . This amounts essentially to finding the invariants 
(symmetry variables) of a given subgroup of the symmetry group admitted 
by Eq.(|l]). A basis set of symmetry variables I{x,y,z,u) for each vector 
field X is obtained by solving the first order partial differential equation 
XI{x,y, z,u) = 0. Below some cases of special physical interest are consid- 
ered. Let us start from the symmetry generator 

X{f) = xd, + ydy - 29,, (10) 

{f = X, f = y). A set of basis invariants is given by 

/ , y 

z = z , t = —, 

X 

W{z, t) = 2\nx + u{x,y,z). (11) 
Then, substitution from (0) into Eq.([l|) yields 

(1 + t^)Wu + 2tWt + 2 = /c(e^),,. (12) 

Exact solutions to Eq. (|1^) can be provided in the framework of group analy- 
sis. In doing so, we have that Eq.(|l2D admits the Lie-point symmetry vector 
fields 

Zo = zd, + 2dw, (13a) 
To = d,, (13b) 
Xq = (1 + t^) arctantc^t — 2(1 + t aTctant)dw, (13c) 
Yo = {l + t^)dt - 2tdw, (13d) 

with [To,Zo] = To, [Yo,Xo] = Y^, 
[To,Xo] = [To,ro] = [Zo,Xo] = [Zo,Yo] = 0. 
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The two-dimensional equation (0) can be reduced to an ordinary differ- 
ential equation in correspondence of any linear combination of the symmetry 
variables (|T3|). For example, the operator generates a configuration which 
can be interpreted as a meron-like excitation. In fact, the reduced equation 
related to is 

(1 + + 2tat = 1 - (14) 

where a = cr(t) = In^ — ^W{z,t), z, t = ^ are invariants. 

Equation (p!4|) can be exactly solved by introducing a function iplt) such 
that 

= 1- Kc-"^". (15) 
Then, Eq.(|l4D can be integrated to give 



{l + t')i,u = 2(c + V^)(l-^i), 



(16) 



where c is an arbitrary constant. 



Now, by setting (j) = t — c — ip, Eq.(|T6]) becomes the Riccati equation 



(17) 



(co= constant) which yields the equation of the hypergeometric type fTO 

Itt + -^-^ It + 77-^ 7 = 
via the transformation 



1+^2 (1+^2)2 
= -(l+t^).-^* 



7 



(18) 
(19) 



The general solution of 



IS 



7 = A 



1 - itX" 



l + it 



+ B 



1 - it 



(20) 



where a = cq/2 and A,B are two constants of integration. 



It 



With the help of (0), (|l|) and (|T8D we have 

4a2 



-2a 



K 



(21) 
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which can be cast into the exphcit form 

AB 

= IG/ta^cos^^ ^, (22) 

in polar coordinates x = r cos 6, y = r sin 6. 

The class of real a corresponds to A = B* and > for k = 1, and to 
A = B* and < for k = —1. Some particular cases are of special interest. 
For example, for a = | and A = B* = 1 + i, from (^) we find (k = 1) 

= i f 1 ± . (23) 

The presence of the sign ± is due to the invariance of Eq.(|lJ) under the 
change of variable t —>■ —t. Consequently, the transformation (p!TD provides 



This expression can be elaborated by resorting to the invariance property of 
Eq.(l) under the conformal transformation (^. Indeed, by setting in ( |24a| ) 
X 2xy, y ^ x"^ — (i.e., by passing to a set of parabolic coordinates), we 
obtain 

= In— , (24b) 

u- = In—. (24c) 
x^ 

Furthermore, performing in u_ the change to parabolic variables, we are led 
to the solution 

«o = ln.^(^ + ^). (24d) 

Now some comments are in order. First, the solution ( |24c| ) corresponds 
to the simplest nontrivial solutions of Eq.([T2|) and (p^, namely = z'^, 
a = (see (|TT])). The solution ( |24bD corresponds to = z^/t'^, a = Int, 
while the configuration (|24d|) corresponds to 

^w^^2l+_^ a = ln^^. (25a) 
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Of course, all the abovementioned expressions for a arise from (^) as par- 
ticular cases. We notice that the related configurations for u{x, y, z) can be 
determined via the conformal invariance property of Eq.(l) starting from the 
simplest solution a = 0. Moreover, by comparing ( P4d|) with ( |24b|) and (|24cj ) 



we have e"+ + e"- 
variable e". 



e"°, which looks "superposition" formula in the 



Second, Eq.(|Tj) plays a role similar to that shown by certain meron equa- 

This can be seen in a heuristic 



tions in the context of Yang-Mills theories [11 



way as follows. For example, let us deal with the solution a 
sociate with it the topological charge defined by [|1^ 

-i g'{t) 

2tT J-oo 



Int, and as- 



(26) 



where g' = and g = e°". Since g' = a'e" , the integral (|26|) is divergent. 
However, it is meaningful as Cauchy principal value, where 1/t is taken as 
a distribution In this sense, q turns out to be zero. Anyway, by lim- 
iting ourselves to consider the integral of g'{t)/g{t) over the interval (— e, e) 
(e > 0), we obtain g = |. This value of q is characteristic of meron-like 
configurations. On the basis of these considerations, we argue that the solu- 
tion M+ = In p-, corresponding to a = Int, behaves as a meron-like excitation 
rather than an instanton-like solution. Furthermore, since the solutions ( |24cj ) 
and ( |24d|) are amenable to u+ via a suitable conformal transformation, 
and Mo can be interpreted as meron-like solutions as well. (Really, through 
the transformations x ^ y, y ^ x, and x 2xy, y x"^ — y'^, from ( |24b| ) 
we obtain ( p4c| ) and (|24d| ) , respectively). 

Another feature confirming the physical character of the configurations 
( |24b| ) , ( p4c| ) and ( p4d|) is their invariance under the transformation x — *■ 
x/(x^ + y^), y —>■ —y/{x'^ + y"^), which corresponds to the inversion rj — *■ 
l/rj. We remind that, as one expects, this property is shared also by the 
gravitational instantons ( ^Sbj ), ( |29a| ) and ( |29bD . Finally, we observe that the 
singularities at x = 0, |/ = of the solutions ( P4b|) , ( |24c| ) and (|24d| ) can be 
shifted through the transformations x ^ x + xq, y ^ y + yo where xq, yo are 
arbitrary constants (see (|^). 

Now, it is instructive to deal with the generators Ji, J2, Pi, P2 of the group 
SL{2,C). 

In correspondence of these operators we obtain four reduced equations of 
the meron-like type, which either coincide formally with ([I^), or are led to 
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(|I2D via a simple change of the independent variable. Hence, these equations 
bring to eight solutions to Eq.(l) which present a meron-like character. For 
brevity, here we shall report only the result for Ji. By using the set of 
invariants 2, = sin y / sinh x , W {z , C,) = m + 2 In sinh x, we find the reduced 
equation {l + C,'^)W^^ + 2C,W^ + 2 = /t(e^)^^, which yields the solutions (k = 1) 



2sinh X 



1 ± 



sm y 
sinh^ X + sin^ y 



(27a) 



Another symmetry generator important by a physical point of view is Y, 
defined by (|5b|). A set of basis invariants are x,y,U{x,y) = u — 2liaz. The 
associated reduced equation is the Liouville equation 



drjdfjU 



K 



whose general solution is 



fviv) fviv) 



(27b) 



(27c) 



for K, = —1, where /(r/) denotes an arbitrary holomorphic function. The 
solution of Eq.(l) corresponding to ( |27c| ) reads 



friiv) fviv) 



1 + ifiv) m 

which produces the Eguchi-Hanson gravitational instanton 

(1 + X2 + 1/2)2' 



(28a) 



(28b) 



by setting /(r^) = 2?7, /(f/) = 2f], {zq = constant). 

From Eq. (|27b| ) other gravitational configurations of physical interest emerge. 
By way of example, if /(?]) = 2c^ri^ , f{jf) = 2c^f/^, we have 



m^z^ 



r J 



N 



N 



(29a) 
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where ro and N are two positive constants such that > 1. The last condi- 



tion ensures the regularity of the r.h.s. of ( 29a ) at the origin and at infinity, 



2N- 



2 and e" 



-2N-2 



respectively. The expression 



where e"_^Q ~ r 

which contains ( |28b|) as a special case, resembles strongly that appearing in 
the context of the self-dual Chern-Simons model developed by Jackiw and 
Pi IjT^. In our framework, the quantity ( |29a|) can be interpreted as an A^- 
gravitational instanton solution where the instantons are superimposed at 
the origin (each with the same scale ro). 

Another nontrivial configuration generated by the symmetry operator 
Y (via the Liouville equation ( p7b| )) is determined by ( |28a| ) for /(r^) = 
2 Y.f=i Ci/ {f] — rfi)^ where q are complex constants. It follows that 



4|E 



N c. 



1 + IE, 



N 



«=1 {ri-rii) I 



2 ' 



(29b) 



which describes separated gravitational instantons depending on AN real 
parameters . 



The metrics related to the solutions (|29a]) and ( |29b|) can be straightfor- 
wardly calculated by Eq.(2). As one expects, the metric associated with ( p9a| ) 
is singularity free and asymptotically flat, while the metric of (|29b|) corre- 



sponds to that devised by Gibbons and Hawking resorting to the multicenter 
ansatz |]15 . 



On the other hand, in opposition to what happens for the gravitational 
instantons, the metric tensors defining the meron-like solutions present sin- 



gularities. For instance, the metric tensor of the configuration (|24c|) is 



{dxf + {dyf +4(rfp) 



-—dy 

X 



dr 



(30) 



with z = p^. The expression (^) has a singularity at x = 0. This is in 
accordance with the meron-like character of ( |24c|) . 

As a final example of remarkable symmetries of Eq.(l), we shall deal with 
the generators i) Vq = Tq + Yq, where Tq and Yq are expressed by ( |13b| ) and 
( |13d| ), and ii) Xq = — h"^ {hdh + Idi — 2du), which comes from (^) 

by setting / = e"'^cosny, f = e"'^ sin ny [h = e"^, / = sinny). 
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The case i) shows the existence of an unsuspected hnk between the heav- 
enly equation (1) and an equation introduced by Schroder in the formulation 
of a combinatorial problem over 100 years ago [0. 

Schroder's equation appears also in the context of the statistical bootstrap 
model and in renormalization theory 0. A set of invariants of Vq is /9 = 
arctant - z/V2, A(/3) = W{f3, ^) - ln(l + (3^). Then, Eq.(pD yields the 
ordinary differential equation 



A/3/3 = K,ie^)/3I3, 

which gives the Schroder equation 

2A - + 1 = /3, 



(31) 



(32) 



for K = 1 and after a suitable choice of the constants of integration. 

The analytical structure of (^2|) was investigated by Hagedorn and Rafel- 
sky [0 who found an integral representation of A in the principal Riemann 
sheet. 

A solution of Eq.(R^) in terms of a power series expansion is given by [H 



(JXJ 

A = Ao- 5]c„(/?o - /?)", 



(33) 



n=l 



where (3o = 21n2 — 1 is the point at which (3{A) has a maximum with value 
Ao = In 2 and the coefficients c„ can be calculated by the recursion relation 



n 



n + 1 



n-1 

CkCn+l~k 

k=2 



(Cl = 1). 



The series converges rapidly in the region < f3 < f3o- Thus, taking account 
of (|33D we get the solution 



u 



In 



n=l 



V2 



(34) 



to Eq.(l), where (3 = 6 — z/\/2, 9 = arctan|, 
W-2\nx. 

The case ii) leads to the reduced equation 



^/x'^ + y"^, and u 



U, 



99 



(35) 
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where z, p = -e sinn?/ and U{q,z) = u + 2nx are invariants. Equation 



(|35D is a two-dimensional version of Eq.(l). Special solutions to Eq. (|35[) are 
U = 21n(2;/cosp) , for k = 1 ITH], and U = 21n (2/coshp) for k = —1, so 
that the functions 

z z 

Ui = 2 In — 2nx, M2 = 2 In — - — — 2nx (36) 

cosp coshp 

solve, in correspondence, the original equation (1). 

At present, the physical interpretation of the configurations and 
( |36D remains to be elucidated, and deserves further investigation. Here we 
remark only that (|3^ ) constitutes, at the best of our knowledge, the only 



known example of a solution to Eq. (1) in which the variable z appears in a 
nontrivial way (usually, the dependence of e" on z looks as a polynomial of 
the second degree). 

To conclude this note, we emphasize the role of the symmetry approach 
in the classification of the gravitational excitations allowed by the heavenly 
equation. 
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